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Abstract 

It is known that Lorentz covariance fixes uniquely the current and the associated guidance law 
in the trajectory interpretation of quantum mechanics for spin-^ particles. In the nonrelativistic 
domain this implies a guidance law for electrons which differs by an additional spin-dependent 
term from the one originally proposed by de Broglie and Bohm. Although the additional term in 
the guidance equation may not be detectable in the quantum measurements derived solely from 
the probability density p, it plays a role in the case of arrival-time measurements. In this paper 
we compute the arrival time distribution and the mean arrival time at a given location, with and 
without the spin contribution, for two problems: 1) a symmetrical Gaussian packet in a uniform 
field and 2) a symmetrical Gaussian packet passing through a ID barrier. Using the Runge-Kutta 
method for integration of the guidance law, Bohmian paths of these problems are also computed. 
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I. INTRODUCTION 



In classical mechanics, each particle follows a definite trajectory, and so it is clear what 
is meant by the time at which a particle arrives at a given place. In quantum mechanics, 
as opposed to classical physics, the meaning of the arrival time of a particle at a given 
location is not evident when the finite extent of the wave function and its spreading becomes 
relevant. Moreover, in the quantum case one expects an arrival time distribution, and there 
are different proposals for this (see e.g. the review article in [1], and the book [2].) As 
pointed out by Hannstein et al [3], these arrival time distributions have been controversial, 
since they are derived from purely theoretical arguments without specifying a measurement 
procedure. The lack of a self-adjoint arrival-time operator conjugate to the free Hamiltonian 
lies at the core of the difficulties in the formulation of quantum arrival times [4]. In Bohm's 
interpretation of nonrelativistic quantum mechanics [5-11], an electron is a particle that has 
a well defined trajectory and so the definition of arrival-time distributions is unambiguous. 
In Bohm's trajectory approach to the calculation of various characteristic times, it is only 
the particle component of the particle that is being clocked. If the arrival-time detector is 
not included in the Hamiltonian, then one has an expression for the ideal or intrinsic arrival- 
time distribution [12]. The arrival-time problem is unambiguously solved in the Bohmian 
mechanics, where for an arbitrary scattering potential V(x), one finds [12-15] for those 
particles that actually reach x = X; the arrival-time distribution being given by the modulus 
of the probability current density, i.e., |J(X,T)|. In nonrelativistic quantum mechanics, the 
form of the current density J is not uniquely determined by the continuity equation, a point 
which has been mentioned by a number of authors [16-18]. It is determined only up to 
a divergenceless vector. For instance, one can construct a new current J' by adding the 
divergenceless current J s to the current J. The newly defined current J' = J + 3 S also 
satisfies the continuity equation, with the same probability density. Although the additional 
contribution in the guidance equation may not be detectable in quantum measurements 
derived solely from the probability density p, it plays a role in the case of arrival-time 
measurements. In the case of spin- 1/2 particles, Holland [19] showed that the particle current 
in the relativistic spin- 1/2 Dirac theory is unique. Demanding that the non-relativistic spin- 
1/2 particle current be obtained from the nonrelativistic limit of the Dirac current, this 
nonrelativistic current is also unique. Struyve et al [20] have considered the uniqueness 
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of paths for spin-0 and spin-1 particles. The spin-dependent Bohm trajectories have been 
investigated for hydrogen eigenstates [21] and an electronic transition in hydrogen [22]. 

It was argued that for free spin eigenstates, spin contributions would in principle be 
experimentally distinguishable for arrival-time distributions of spinless and spin-1/2 particles 
[23] . Holland et al [24] have explored some of the implications of the modified guidance law 
in the case of the two-slit quantum interference. In this case, the trajectories cross each 
other and they also cross the symmetry axis. 

The aim of the present paper is to explore some of the implications of the revised nonrel- 
ativistic guidance equation. Sec. II contains a very brief review of relevant parts of Bohm's 
interpretation of quantum mechanics. In Sec. Ill, by using the nonrelativistic limit of the 
Dirac current density, numerical computations of the effect of the spin-dependent term on 
arrival time at a given location, are presented for both a Gaussian wave in a uniform field 
and a Gaussian wave passing through a barrier. 

II. BOHM'S TRAJECTORY INTERPRETATION OF QUANTUM MECHANICS 
AND THE ARRIVAL TIME DISTRIBUTIONS 

In nonrelativistic Bohmian mechanics the world is described by point-like particles which 
follow trajectories determined by a law of motion. The evolution of the positions of these 
particles are guided by a wave function which itself evolves according to the Schrodinger 
equation. Bohmian mechanics makes the same predictions as the ordinary nonrelativistic 
quantum mechanics for the results of any experiment, provided we assume a random distri- 
bution for configuration of the system and the apparatus at the beginning of the experiment, 
given by p(x, 0) = ^(x, 0)\l/(x, 0). If the probability density for the configuration satisfies 
p(x, to) = ^(x, £ )\I/(x, t ) at some initial time to, then the density to which this is carried 
by the continuity equation at any time t is also given by p(x, t) = \E^(x, t)\I/(x, t) [25]. 
As most of the quantum measurements boil down to position measurements, Bohm's the- 
ory and the standard quantum mechanics generally yield the same detection probabilities. 
The situation is different, however, if one considers, for example, measurements involving 
time-related quantities, such as arrival times, tunnelling times, etc. Bohm's theory makes 
unambiguous predictions for such measurements, but there is no consensus about what these 
quantities should be in the conventional quantum mechanics [26] . Given the initial position 
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x(°) = x(t = 0) of a particle with the initial wave function \&(x, t = 0), its subsequent trajec- 
tory x(x'°',t) is uniquely determined by the simultaneous integration of the time dependent 
Schrodinger equation, and the guidance equation = v(x(£),i), in which v = -. For 
spinless particles, the Schrodinger equation yields: 



J(x,t) = (h/m)Im(il>*(x,t)Vil>(x,t)) . (1) 

But, for systems with more than one spatial dimension the probability current density, and 
hence the particle equation of motion, is not uniquely defined within nonrelativistic quantum 
mechanics. But, Holland [19] has shown that the probability current density deduced from 
the continuity equation is uniquely defined for Dirac electrons, when Lorentz covariance is 
imposed. Taking the nonrelativistic limit for a spin eigenstate in the absence of a magnetic 
field, one gets 



J(x,t;s) = J(x,t) + (fi/2m)(Vp(x,t))xs 

= (h/m)[Im(ip*(x,t)Vip(x,t)) +ite(^*(x,t)VV>(x,t)) x s] , (2) 



where 



-s = -xVx , (3) 



is the spin vector associated with the spin eigenstate x> which is a two-component spinor 
normalized to unity (x^X = !)• The nonrelativistic particle is, in effect, guided by the 
two-component wave function \P(x, t;s) = ip(x,t)x, and the original expression for the 
nonrelativistic velocity field should then be replaced by 



v(x,t;s) = ^§. (4) 
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A striking property of the spin-dependent term is that the components of the particle 
motion in orthogonal directions are generally mutually dependent, even when the wave 
function factorizes in these direction. For a 3D system with a wave function which is in the 
factorized form t/>(x, 0) = i/) x (x, 0)ip y (y, 0)if) z (z, 0) at t — 0, we have 

V>(x,t) = il) x (x,t)ipy(y,t)il) e (z,t), 

at t > 0, and for 's = (0, 0, 1), Eq. (4) becomes 



_ [/m( _ ) + Jfe( _ 

ft rr^/^llN n./da^I 



% = - [j m (-«p)-ite v 

«, = - , (5) 

where we have used the notation d x ip x = and so on. For a nonrelativistic guiding wave 
of the form \&(x, t;"s) = ip(pc,t)x, with ip('x,t) being a solution of the Schrodinger equation 
and x being a fixed spinor, Bohm's trajectory result for the distribution of particle arrival 
times T, at a given location x = X, is given by [12] 



POO 

n(T,X;s) = |J(X,T;s)|/ / d*|J(X,*;s)| . (6) 

Jo 

The mean arrival time of the particles reaching a detector located at X is given by 



dt t |J(X,t;s)| 

(7) 



/ dt t II(t,X;s) = 
Jo I 

Jo 



dt |J(X,t:s)| 



in the presence of the spin-dependent contribution and 



5 



f 

Jo 

J 

Jo 



oc 



I 



oc 



dt t \J(X,t)\ 



dt t Ui(t,X) 



(8) 



dt |J(X,*)| 



in the absence of the spin-dependent contribution. It must be mentioned that the definition 
of the mean arrival time used in Eq. (8) is not a uniquely derivable result within standard 
quantum mechanics. We compute the arrival time distribution and the mean arrival time 
at a given location for two problems: 1) a symmetrical Gaussian packet in a uniform field 
and 2) a symmetrical Gaussian packet passing through a ID barrier with and without the 
spin contribution. Using the Runge-Kutta method for the integration of the guidance law, 
Bohmian paths of these problems, have been also computed. 

A. Symmetric Gaussian packet in a uniform field 

For the potential V(x) = K.x, the wave function and the probability density at time t, 
calculated from [8], are given by: 



provided that this symmetrical Gaussian wave function is centred around the origin at t — 0. 
In the Eq. (9), s t = a (l + iht/2mal), a = a [l + (^r) 2 ] 1 ^ 2 and u is the group velocity. 
The Eq. (9) is in the factorized form in three coordinate directions. We choose the uniform 
field and the group velocity to be in the x direction i.e., ^(x) = K.x = Kx and u = (u, 0, 0) 
and the spin vector in the z direction. The speed of the wave packet's centre in the y and 
z directions has been set equal to zero and therefore this function will spread but will not 
propagate in these directions. Using the Eq. (5), one gets for this factorized wave 



V>(x, t) = (27rs t 2 )- 3/4 exp{-(x - ut + Kt 2 /2m) 2 /4s t a 

+ (im/h)[(u - Kt/m).(x - -ut) - K 2 t 3 /6m 2 ]} , 
p(x,t) = (27ra 2 )~ 3/2 exj9-(x-ut + Kt 2 /2m) 2 /2a 2 , 



(9) 
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Kt h 2 t Kt 2 h 

h 2 t h Kt 2 

h 2 t 

Vz = 4m*a* + hH2 Z ' (10) 

from which one can compute Bohmian trajectories. In the z direction, one gets z(t) = 
zq(1 + 4^~? ) 1 ^ 2 ) in which is the ^-component of the initial position of the particle. 



B. Symmetric Gaussian packet passing through a barrier in the x-direction 

suppose V(x) = V(x) = Vq for < x < d and zero otherwise. The wave function at time 
t is given by 



i/} x (x,t) = / dk(j) k (x)ijj(k)e 



-iE k t/h 



*l> v {y,t) 



1 



(27^)1/4 

1 

WW 1 



e 



y 2 /is t a 



1>'M = TtAkUI^ 2 /^ , (11) 



in which 



k(x) = Mx) = -^T{k)e lkx x>d 
\/2tt 



$( k ) = ( ^0)l/4 e - CT g( fc - fc o) 2 ) (12) 



where u = (hko/m, 0, 0) is the group velocity and 

Tk _ e -ikd ' 2 M h 2 k 2 /2m ^ V , 

2kqcos(qd) — i(q 2 + k 2 )sin(qd) 

T k = — A— e"** h 2 k 2 /2m = V , (13) 

2 + ikd 
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is the transmission probability with q = 2m (H 2 k 2 / 2m — Vq) / h. Again, the speed of the 



wave packet's centre in y and z directions has been set equal to zero and therefore this 
function will spread but not propagate in these directions. For the motion in the z direction, 
one gets 



in which z is the z-component of the initial position of the particle. 

III. NUMERICAL RESULTS 

Numerical calculations are presented for a case of symmetric Gaussian packet in a uniform 
field and a case of symmetric Gaussian packet passing through a barrier in the x-direction. 
We compute the mean arrival time and the arrival time distribution with and without 
the spin-dependent term. Using the Runge-Kutta method for the integration of guidance 
equation, we compute some Bohmian paths for each case. 

A. Symmetric Gaussian packet in a uniform field 

The initial wave packet is peaked at (xq = 0, yo — 0, z = 0), with o~o = 5A. The detector 
position is chosen at (x = 20, y = 20, z = 20)A. 

Fig. 1 shows the arrival time distribution of electrons for K = mg, corresponding to 
particles in a gravity field, with Eq = \mu 2 = 5ev and g = 9.8m/ s 2 . 

Fig. 2 shows the mean arrival time versus the mass of the arriving particle for a fixed 
value of Eq (Eq = 5ev). From this figure, it follows that the spin-dependent term increases 
the mean arrival time (r > r,) and has very small effect on the mean arrival time at the 
given location, and the quantity r — n is decreased with mass. 

Fig. 3 shows the mean arrival time of electrons versus the ratio of the group velocity to 
the light velocity u/c (Note that we have fixed the mass). It follows from this figure that 
the spin-dependent term increases the arrival time (r > Tj) and has very small effect on the 
mean arrival time at the given location and that the quantity r — is increased with the 




(14) 
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initial group velocity at first and then is decreased with it. 

Fig. 4 shows some Bohmian paths for electrons. One can see that some trajectories cross 
each other when one considers the spin-dependent term. It should be noted that trajectories 
cross each other for different values of the y coordinate. Consequently there is no problem 
with the single-valuedness of the wave function. The spin-dependent term can also change 
the fate of the individual trajectories. There are some trajectories that don't reach the 
detector location, with the spin-dependent term, but reach the detector location without 
this term. 

B. Symmetric Gaussian packet passing through a barrier in the x-direction 

The initial wave packet is peaked at {xq = — 10o"o, £/o = 0, zq = 0) with o"o = 5A. The 
detector location is chosen at (x = 20, y = 20, z = 20)A. 

Fig. 5 shows the arrival time distribution of electrons, for Vq = 8ev, Eq = ^mu 2 = lOev 
and d = lOA. 

Fig. 6 shows the mean arrival time of electrons at the given location, versus width of the 
barrier. From this figure, it follows that the spin-dependent term increases the mean arrival 
time (r > Tj) and the quantity r — Tj is increased with the width of the barrier. 

Fig. 7 shows the mean arrival time of electrons versus the ratio of the group velocity 
to the light velocity u/c for d = 10A. It follows from this figure that the spin-dependent 
term increases the arrival time (r > Tj) and has very small effect on the mean arrival time 
at the specified location, and that the quantity t — Tj is decreased with the initial group 
velocity. Note that we could not decrease the initial energy (or the initial group velocity) 
as we want, as in the Fig. 3, because if Eq = \mu 2 <C V , then the transmission probability 
would be very small and no particle can cross the barrier to reach the detector position. 

Fig. 8 shows some Bohmian paths for electrons. One can see that trajectories cross 
each other (Fig. 8), when one considers the spin-dependent term. As the Fig. 4 they 
cross each other for different values of the y coordinate. The spin-dependent term can also 
change the fate of the individual trajectories. There are some trajectories that don't reach 
the barrier, with the spin- dependent term, but they cross the barrier without this term. 
However, ensemble averaging yields the same transmission and reflection probabilities. The 
spin-dependent term plays a role only in time related measurements. 
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IV. SUMMARY AND DISCUSSION 



Spin has a very small effect on the mean arrival time at a given location in the non- 
relativistic domain, but it has a significant effect on the fate of individual Bohmian paths. 
Bohmian paths can cross each other, when one considers the spin-dependent term. The 
quantity r — r, changes with the mass of the arriving particle, with the group velocity of the 
wave packet and also with the width of the barrier, but, the difference is very small. 
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Figure captions 

Figure 1: The arrival time distribution of electrons versus time, without the spin-dependent 
term (solid curve) and with the spin-dependent term (shown by triangles), for a symmetric 
Gaussian packet in a uniform field. 

Figure 2: a) The mean arrival time versus the ratio of the mass of the arriving particles to the 
mass of electron, without the spin-dependent term (solid curve) and with the spin-dependent 
term (shown by triangles) and b) the difference between T{ and r versus m/m e for a symmetric 
Gaussian packet in a uniform field. 

Figure 3: a) The mean arrival time of electrons versus the ratio of the group velocity to the 
light velocity, without the spin-dependent term (solid curve) and with the spin-dependent 
term (shown by triangles) and b) the difference between Tj and r versus u/c for a symmetric 
Gaussian packet in a uniform field. 

Figure 4: Some Bohmian trajectories of electrons in x-t plane with the initial position 
(x^ , = yj 2<7q — x(°) 2 , = 0) a) in the absence of the spin-dependent term and b) 
in the presence of the spin-dependent term, for a symmetric Gaussian packet in a uniform 
field. 

Figure 5: The arrival time distribution of electrons versus time, without the spin-dependent 
term (solid curve) and with the spin-dependent term (dashed curve), for a symmetric Gaussian 
packet passing through a barrier in the x-direction. 

Figure 6: The mean arrival time of electrons versus the width of the barrier, without the 
spin-dependent term (solid curve) and with the spin-dependent term (dashed curve), for a 
symmetric Gaussian packet passing through a barrier in the x-direction. 
Figure 7: The mean arrival time of electrons versus the ratio of the group velocity to the light 
velocity without the spin-dependent term (solid curve) and with the spin-dependent term 
(dashed curve), for a symmetric Gaussian packet passing through a barrier in the x-direction. 
Figure 8: Some Bohmian trajectories of electrons in x-t plane with the initial position 
( x (o) ^ y(o) _ y/ 50<Tq — x(°) 2 , z^ = 0) a) in the absence of the spin-dependent term and b) 
in the presence of the spin-dependent term, for a symmetric Gaussian packet passing through 
a barrier in the x-direction. 
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FIG. 1: The arrival time distribution of electrons versus time, without the spin-dependent term 
(solid curve) and with the spin-dependent term (shown by triangles), for a symmetric Gaussian 
packet in a uniform field. 
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FIG. 2: a) The mean arrival time versus the ratio of the mass of the arriving particles to the 
mass of electron, without the spin-dependent term (solid curve) and with the spin-dependent term 
(shown by triangles) and b) the difference between r, and r versus m/m e for a symmetric Gaussian 
packet in a uniform field. 
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FIG. 3: a) The mean arrival time of electrons versus the ratio of the group velocity to the light 
velocity, without the spin-dependent term (solid curve) and with the spin-dependent term (shown 
by triangles) and b) the difference between n and r versus u/c for a symmetric Gaussian packet 
in a uniform field. 
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FIG. 4: Some Bohmian trajectories of electrons in x-t plane with the initial position (x^°\y^ = 

y^Ug — x(°) 2 , = 0) a) in the absence of the spin-dependent term and b) in the presence of the 
spin-dependent term, for a symmetric Gaussian packet in a uniform field. 
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FIG. 5: The arrival time distribution of electrons versus time, without the spin-dependent term 
(solid curve) and with the spin-dependent term (dashed curve), for a symmetric Gaussian packet 
passing through a barrier in the x-direction. 
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FIG. 6: The mean arrival time of electrons versus the width of the barrier, without the spin- 
dependent term (solid curve) and with the spin-dependent term (dashed curve), for a symmetric 
Gaussian packet passing through a barrier in the x-direction. 
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FIG. 7: The mean arrival time of electrons versus the ratio of the group velocity to the light 
velocity, without the spin-dependent term (solid curve) and with the spin-dependent term (dashed 
curve), for a symmetric Gaussian packet passing through a barrier in the x-direction. 
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FIG. 8: Some Bohmian trajectories of electrons in x-t plane with the initial position (x^°\y^ = 

y^OcjQ — x(°) 2 , = 0) a) in the absence of the spin-dependent term and b) in the presence of the 
spin-dependent term, for a symmetric Gaussian packet passing through a barrier in the x-direction. 
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